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ABSTRACT 

Langlands et al. considered two crossing probabilities, ir^ and tt^v, in their ex- 
tensive numerical investigations of critical percolation in two dimensions. Cardy 
was able to find the exact form of tt^ by treating it as a correlation function of 
boundary operators in the Q — ► 1 limit of the Q state Potts model. We extend 
his results to find an analogous formula for irhv which compares very well with the 
numerical results. 
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1 Introduction 



Critical percolation in two dimensions has been intensively studied numerically by Langlands 
et al., |l|, One of the results they find is that various crossing probabilities, tt^, 717^ (to 
be defined below) are invariant under conformal transformations of the plane. In ||, by 
combining a standard identification of critical percolation and the 1 state Potts model with 
boundary conformal field theory techniques, Cardy was able to find the exact form of 717^ 
There are three steps needed to extend this to 717^: 1. Show that 717^ corresponds to some 
particular boundary conditions in the 1 state Potts model, and thus is given by a correlation 
function of conformal field theory boundary operators; 2. identify theses boundary operators; 
3. deduce a differential equation for the correlation function, and compare the result with 
■Khv- in this letter we complete steps 1 and 3. Of course this is far from being a complete 
derivation of 717^, but having shown that it is given by a correlation function of boundary 
operators, then one knows that 717^ must satisfy one of a discrete set of known differential 
equations, and examining the simplest of these we find a single function which has suitable 
asymptotic behaviour. Comparing this function with the numerical data for 717^ presented in 
]|], our agreement is as good as that of Cardy's formula for 717^ 



2 Critical percolation and boundary conformal field theory 



For an extremely good presentation of critical percolation, and the full content of the hy- 
potheses of universality and conformal invariance, we refer the reader to [jij. At its simplest, 
one can envisage a rectangle, and take as observables all combinations of probabilities of 
crossing from one set of points on the boundary to another disjoint set. The hypothesis of 
conformal invariance is that such probabilities are invariant under conformal transformations 
of the plane. For example, the probabilities, 717^ of crossing from the left side to the right, and 
Tihv of being able to cross simultaneously from left to right and from top to bottom, depend 
only on the aspect ratio r of the rectangle. It is easy to see that, with r = height/width, 

Tr hv (r) = ir hv (l/r) , n hv (r)/Tr h (r) -> 1 as r -> . (2.1) 

Langlands et al. provide numerical values for 7^ and 717^ for aspect ratios 0.136 < r < 7.351. 

The property of the crossing probabilities described above, that they are invariant under 
conformal mappings of the plane, is that of correlation functions of h = boundary operators 
in a c = conformal field theory. (For a review of boundary operators in conformal field 
theory, see e.g. [Q, ||]). It was in this way that Cardy was able to find a formula for 717,, (r). It 
will be instructive to repeat the relevant points of his derivation here. 
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The combination (1 — TTh) may be identified with the Q — > 1 limit of the Q state Potts 
model partition function with specific boundary conditions: namely free boundary conditions 
on the vertical sides of the rectangle, and the spins fixed to different values on the horizontal 
sides. One can identify this partition function as the correlation function of four operators, 
placed at the corners of the rectangle, which change the boundary conditions. 

By considering the 2 and 3 state Potts models, for which the correspondence with c = 1/2 
and c = 4/5 conformal field theories is very well known, Cardy could identify the relevant 
boundary changing operators as type (1,2) boundary primary fields. This implied that the 
correlation function of four such operators should satisfy a second order differential equation. 
This equation is very easy to solve, and by taking the appropriate combination of the two 
independent solutions of this equation, Cardy could fit the data of [P to a very high degree 
of accuracy. In fact, this fit was so good, that Langlands et al. decided to adopt Cardy's 
result as TTh, rather than their own numerical data, in their later paper on more complicated 
calculations Q. 

3 Boundary conditions for n^v 

One can similarly find TThv from the Q — > 1 limit of the Q state Potts model partition function 
with certain other boundary conditions which, exactly as for n^, are not themselves defined 
for Q = 1. There is a standard argument relating the partition function in the Potts model 
to a sum over a graphical expansion Q. In this, one considers all sets of 'clusters' £, where 
a cluster is a connected set of bonds on a lattice, and no two clusters intersect. Then the 
partition function may be written as a sum over clusters, with weights which depend on the 
type of cluster, corresponding to the type of boundary conditions one wishes to impose. For 
our purposes, it will be enough to consider a partition function defined for a square lattice 
on a rectangle with the sides of the rectangle consecutively numbered 1 to 4, and to consider 
an expansion of the form 

Z(a,b,C,d,Q) = J2p m (± -p) B - B (O a Na(O b N^) c N^) d N-N a (^N b (0-N^) ? (31) 

where p is the normalised Boltzmann weight for an edge connecting sites with the same spin, 
and which will be set to the critical probability p c = 1/2 for a bond to be open, when we finally 
take Q — > 1. In ( p.l[ ), B is the total number of bonds, -B(£) is the number of bonds in £, N is 
the total number of clusters, N a , iVj, and N c are the numbers of clusters intersecting no sides of 
the rectangle, exactly 1 side and two adjacent sides respectively, and a, b, c and d are weights 
which depend on the boundary conditions. Taking a = b = c = d = Q one recovers the 
Q state Potts model partition function for free boundary conditions, and Z(l, 1, 1, 1,1) = 1 
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in the 1 state Potts model. However, if one can arrange boundary conditions such that, as 
Q — ► 1, that a — ► 1, b — ► 1 and c — ► 1 and d — > 0, then one will recover the probability that 
no cluster intersects non-adjacent boundaries, i.e. 1 — 7r/j — ir v + tt/h;. 

Let us consider the following boundary condition, which we denote by {AB . . . C): in this 
boundary condition the spin at each site in the boundary may take its value freely in any of 
the states A, B, . . . C. If we assign the conditions bi to the four sides of the rectangle in turn, 

b\ = {A, B, Xn,X\2, ■ ■ ■ , X\ n ) , 62 = (B, C, X21, X22, ■ ■ ■ , X2n) , 

63 = (C, D,Xs\, X32, ■ ■ ■ >^3n) , b^ = (D, A,X^i, X42, ■ ■ ■ , X^n) , (3.2) 
where all the spins {Xij} and {A, B, C, D} are distinct, then we see that the partition function 



of the Q-state Potts model with these boundary conditions will be in the form (3.1) with 
weights a = Q, b = n + 2, c=l and d = 0. Then, if we take n = (Q — 5)/4 (which is perfectly 
valid for Q = 5, 9, . . .) and then Q — > 1, we will recover the required weights for our crossing 
probability, i.e. 

hin Z(Q, (Q + 3)/4, 1, 0, Q) = 1 - Tr h - n v + n hv = ir hv , (3.3) 

where we take for granted the universality of critical percolation and set 7r^ + ir v = 1 as is the 
case for percolation by sites on a triangular lattice. 

Unfortunately, the boundary changing operator which changes the boundary conditions 
from type (A, B, X\\ , X\2, ■ ■ ■ , X\ n ) to (B, C, X21 , X22, • • • , X2 n ) does not exist in the 2 state 
Potts model for any n, and has only been identified to date for n = in the 3 state Potts 
model [||, and consequently we can not deduce a differential equation for Z. However, on 
the basis of this argument, we can certainly expect that irh v is expressible as some correlation 
function of four h = conformal boundary primary fields. 



4 Differential equations for correlation functions 



Correlations of conformal primary fields (fth( z ) can satisfy differential equations, one differen- 
tial equation for each null state in the highest weight representations of the Virasoro algebra 
with highest weight state 0^(0) jO) , and the general method by which such differential equa- 
tions are derived is given in M. Since the space of null vectors is a highest weight space, one 
can restrict attention to the differential equations arising from the vanishing of highest weight 
states. For h = c = it is easy to compute the highest weight states with low conformal 
weight in the Verma module and we give the first four of these in table |l]. The four-point 
functions of four h = boundary primary fields will only depend on the cross ratio r), of their 
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coordinates, so we will write 



(<t>(z 1 )<j>(z 2 )<t>(z3)<j>(z4))=F(r ] ) 

where r\ = ((z^ — z/±){z2 — zi)/((z3 — Z\){zi — Z4)). The differential equations arising from the 
first three null vectors are equally easy to derive, and are given, along with their solutions, in 
table |2[ The conditions ( |2.1| ) translate into 

F(n)=F{i-n), 3 Fn ^/o A/ o ^ WT7d as ^°- ^ 

t] ' a-Fi(l/3, 2/3; 4/3; 77) r(l/3) 2 

It is impossible to find a solution of type or F^ which satisfies but there is a 

unique solution of type F^ which does, and which can be presented variously as 

r(2/3) f* dt 2 f z dt [* du 

J~\Z) 



r(l/3) 2 7 - i)]2/3 3T(l/3)r(2/3) Jo -t)]2/3 7 [ u (i-„)]i/3 
r(2/3) r " 1 r <" -2^1(1, 4/3; 5/3; t) 



r(i/3) 2 7 - 1)]2/3 r(i/3)r(2/3) 7 - i)] 2/3 

= ^0^ 1/3 2^(1/3, 2/3; 4/3; r(1/3) r (2/3) 3^(1, 1, 4/3; 2, 5/3; z), (4.2) 

where 2^1 is the standard hyper geometric function, and 3-F2 is a generalised hypergeometric 
function. 



5 Comparison with data and discussion 

The case treated in |2| is of the crossing probabilities defined for a rectangle of aspect ratio 
r. If we take the points Zj to be at (—A; -1 , —1, 1, A; -1 ), then we can map these to the corners 
of a rectangle by 

f z dx 

ID = I — — 

Jo ^{l-x 2 ){l-k 2 x 2 ) ' 
in which case the aspect ratio of the rectangle is given by r = K(l — k 2 )/2K(k ), where 
K (u) is the complete elliptic integral of the first kind. The cross ratio of these four points is 
77 = ((1 - fc)/(l + k)f. 

In figure 1 we plot log(7T/j t) ) against log(r) for both the numerical results obtained in E| 
and for the function J-{x) above. The agreement is excellent^]. 

There are still several points which remain unclear in the results presented here. Firstly, 

the integral expression of ir^ found by Cardy can be easily identified as one of the standard 

solutions of Dotsenko and Fateev Q. However, despite the fact that T satisfies the standard 

1 In fact Langlands et al. only computed -k^v for r < 1, since it is invariant under r — > 1/r, and consequently 
half of the data points we give are duplicated. 
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differential equations for a type (3, 4) primary field, we have not yet been able to find a 
way to express ( |4.2| ) in Dotsenko-Fateev form. If one allows more general vertex operator 
constructions, such as for VK-algebras, then it is easy to see how the troublesome second term 
in T might arise as a screened four-point function: one way, using purely bosonic vertex 
operators, is 

lim [*dt fdu (3a + P + 7 | V a (w) V a {l) V (z) V p {t) VJu) V a (0) |0) , 

w->oo Jo J 

where a 2 = 1/2, a-(3 = —2/3, 0-7 = —1/3, = and where V\(z) =: exp(i\ ■ 4>{z)) : 
is a bosonic vertex operator, V a represents a weight field, and J Vp and / F 7 are screening 
charges. However, we have not yet found any pressing case for such an interpretation for any 
particular extended algebra. 

Secondly, although T{z) satisfies the fifth order differential equation in table ^, it also 
satisfies the third order differential equation 

± (z(z-l))V3 A. ( z ( z -l)f/3 * jr =0 . (5.1) 

az az az 

This is a much more appealing equation as the three independent crossing probabilities 



l,TTh,7Thv span the solutions to (5.1), but this equation does not arise from the vanishing 
of a vector at level 3 in the h = Virasoro Verma module. It would be very nice indeed 
if this equation could be derived from a null-vector vanishing for some generalisation of the 



Virasoro algebra, such as those in [hd, 11], but we have not been able to do this yet. 



6 Conclusions 



From a simple extension of Cardy's method of Q, we have found an excellent candidate 
for the crossing probability irhv which was investigated numerically by Langlands et al. in 
0. Unfortunately, we have not been able to find a real derivation of our result, but work 
is in progress which should elucidate rather more the nature of the conformal field theory 
underlying this result. 
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Table 1: Low lying highest weight vectors in the h = c = Verma module. 



level n 


Null vector M{n) 


1 

2 
5 

7 


1 

( 
( 


-i|0) 

L_!L_! - (2/3)L_ 2 ) |0> 

L 3 _ x - 6L_ 2 L_i + 6L_ 3 ) JV(2) 

^ Llj - (40/3)L_ 2 L 3 _ 1 + (256/9)L_ 2 L_ 2 L_i + (52/3)L_ 3 L_iL_i ^ 
^ -(256/9)L_ 3 L_ 2 - (104/3)L_ 4 L-i + (208/9)L_ 5 y 


AT(2) 



Table 2: Differential operators and solutions 



level n 




1 

F m 


A 

dz 
Cl 


2 

^(2) 


(*(*-l))- 2/3 £ (^-1)) 2/3 

Cl + C2 7o (x(l-x)) 2 /3 


A 

dz 


5 
i?(5) 


(z(.-l))- 2 £ (z(z- 
f z dx 


l))4/3 
C 2+ J 


£ (z(z - I)) 2 / 3 £ 
f x c 3 + c A t + c 5 t 2 ' 




Cl + Jo (x(l - X)) 2 / 3 


(t(l-t))^ _ 
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